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ABSTRACT 

Accurate measurement of the cluster mass function is a crucial element in efforts to constrain 
structure formation models, the normalisation of the matter power spectrum and the cosmo- 
logical matter density, and the nature and evolution of dark energy. Large weak lensing surveys 
of ~ 20, 000 galaxy clusters and groups will be key tools in the observational pursuit of that 
goal; first-generation surveys such as those using the Sloan Digital Sky Survey are already be- 
ginning to calibrate previous measurements of the mass function from X-ray observations and 
to extend existing constraints down to galaxy group scales. These weak lensing studies pro- 
ceed by stacking the lensing signals of many clusters and groups binned by mass-correlated 
observables such as richness and luminosity; typically such analyses ignore the triaxial struc- 
ture of dark matter halos on the assumption that the averaging of many shear signals within 
each mass bin makes its effects (as large as factors of two in mass model parameter estimates 
from individual clusters) negligible. We test this assumption utilizing a population of 15,000 
analytic triaxial dark matter halos spanning two dex in mass, and find that triaxiality can bias 
3D virial mass estimates compared to those for a spherical population by a few percent if 
suboptimal mass estimators are used. This bias affects not only direct lensing constraints on 
the mass function but can also affect the scatter and normalization of the mass-observable re- 
lations derived from lensing that are so crucial to constraining the cluster mass function with 
large samples. However, we demonstrate that a careful choice of mass estimator can remove 
the bias very effectively if the lensing signals from a sufficient number of triaxial halos are 
averaged together, and further quantify that sufficient number for adequate shape averaging. 
We thus show that by choosing observable bins to contain an adequate number of halos and 
by utilizing a carefully chosen 3D mass estimator stacked weak-lensing analyses can give 
unbiased constraints on the triaxial mass function. 

Key words: gravitational lensing - cosmology: theory - dark matter - large-scale structure of 
universe - galaxiesxlusters: general - methods: statistical 



1 INTRODUCTION 

The number of galaxy clusters and groups present in the uni- 
verse as a function of mass and redshift is a strong constraint 
on cosmological models: measuring this cluster mass function 
over a wide range of masses and redshifts is a crucial element 
in the effort to constrain structure formation models, cosmolog- 
ical parameters Qm and erg, and the nature and evolution of 
dark energ y. The mass f unctio n has been well studied in simula- 
tions (e.g. Evrard et alj J2002f)) and in nature i n the X-ray (e.g. 
Ivikhlinin et alj d2008h ; jReiprich & Boehringerl J2002h ). the opti- 
cal (e.g. iBahcall et al. I <t2003|)), and in first-generation weak lens- 
ing stu dies in large surveys such as the Sloan Di gital Sky Survey 
(SPSS, lYork etafl feOOCb ; Ijohnston et al] [2008)). A detailed re- 
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view of the theory of the cluste r mass function and its observational 
pursuit is given by Voit (2005). 

Weak lensing is of growing importance in the effort to con- 
strain the mass function as the next generation of weak lensing sur- 
veys promise much larger samples of galaxy clusters and groups 
over a wider range of redshifts than have been studied to date. It 
may be a preferred method because the mass measurements it pro- 
vides do not depend on baryonic physics and are largely indepen- 
dent of the observables by which lens systems are sorted into mass 
bins (X-ray temperature or luminosity; optical luminosity or rich- 
ness), making it less likely to be biased by s ystematic errors specific 
to one type of observational data (see e.g. [Mandelbaum & Seljak 
(2007)). However, lensing is subject to its own systematic and sta- 
tistical errors, which must be controlled for it to become reliable 
as a measure of the mass function. The Dark Energy Task Force 
(DETF) sets the goal that the RMS error in the mean and variance 
of the mass estimates in current and next-generation weak lens- 
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Figure 1. Offsets in masses recovered, averaged across 50 noise realizations, for lensing by identical populations of 15,000 spherical and triaxial lensing halos, 
binned using SDSS observables. The top left panel shows the percent offset of the mean recovered lensing mass compared to the true mean mass (Af) 2 oo 
for each luminosity L200 bin. The top right panel is the same but compares with the true median mass med(M2oo)- The bottom panels show the same for 
the richness N200 bins. The solid black lines plot the results for lensing by a population of spherical halos, the blue dot-dash lines show them for lensing 
by the corresponding triaxial halo population, and the red dashed line shows the results for lensing by the triaxial population compared to the true mean and 
median effective spherical mass Mes m eacn Dm - The effective spherical mass Mes is the triaxial mass measure most closely recovered by a stacked lensing 
analysis. The error bars give the standard deviation across the 50 realizations. 



ing surveys be reduced to n o more than 14% and 11% respectively, 
and ideally to less than 2% I Albr echt et alj j2006h) in order to accu- 
rately construct a mass function. A full understanding of all sources 
of error is crucial to achieve those levels. Some possible sources are 
unaccounted-for cl uster substructur e , scatter from lar ge scale line- 
of-sight structures JPodelsorJ d2004l) ; lHoekstr j d2003h >. incomplete 
cluster and group samples, and mischaracterization of the scatter 
and normalization in the mass-observable relations used to bin and 
stack the weak lensing signals. 

In this paper we address one potentially important systematic 
error: the neglect of the triaxial structure of the dark matter ha- 
los that shape clusters an d groups. CDM str uctur e formation simu- 
lations such as those of lShawetall d2006h and iBett et all d2007l) 
indicate populations of significantly triaxial halos that are more 
prolate than oblate with minor axis ratios as small as 0.3. Sev- 



eral studies have examined the importance of this triaxial struc- 
ture in lensing analyses of individual clusters, which typic a lly em - 
ploy spherical, or at best elliptical, models dOguri et al.l (f2005h: 
ICorless & King] d2007h : ICorless & Kind i l2008h : iMeneghetti et al.1 
and found its neglect to cause errors of up to factors of 
two in estimates of virial mass and halo concentration. 

Typically triaxiality is neglected in stacked weak lensing anal- 
yses measuring halo masses because it is argued that the aver- 
aging of the lensing shear signals of many halos randomly ori- 
ented in each bin makes the spherical mass derived from the mean 
signal an accu r ate es timator of the true bin mean mass (see e.g. 
Ijohnston et ail fc(X)8l» . However, Corless & King (2007; 2008) 
showed that while indeed the mean mass across a population of 
triaxial halos fit individually with spherical models is a good esti- 
mator of the mean population value, the scatter is asymmetric and 
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Figure 2. The mean cumulative mass function for halos of M > 1.0 X 
lO 13 Af0 recovered from the 50 noise realizations of Figure [T] for lens- 
ing by identical populations of 15,000 spherical and triaxial lensing halos, 
binned using the SDSS richness observable N200 ■ At low masses the im- 
pact of the triaxiality bias is very small; at high masses, the region richest 
in information about cosmological parameters, it is small but potentially 
significant. 



the error estimates incorrect. Given this, the interpretation of the 
mean shear signal in a bin stretched across a steeply sloped mass 
function is complex and not obvious. Neglecting this complexity 
in mass modelling may lead to mass estimates that cannot be in- 
terpreted as theoretical quantities such as the virial mass, an inter- 
pretation which is necessary to make comparisons between obser- 
vational results and theoretical predictions. The importance of the 
mass function for cosmology and the increasing efforts that will be 
devoted to constraining it in future lens surveys make it crucially 
important to fully understand any impact triaxiality may have on 
the accurate and clear measurement of galaxy cluster and group 
masses. 

To this end, in this paper we directly model the impacts of tri- 
axiality on the masses recovered from weak lensing surveys, using 
an analytic lens population based on an SDSS sample. Section|2]de- 
scribes our simulations of lensing by a large population of triaxial 
galaxy clusters and groups, Section|3]outlines the mass-observable 
binning and simple stacked weak lensing techniques we employ to 
fit their masses, and Section [4] presents our results and some dis- 
cussion. 



2 SIMULATIONS: LENSING THROUGH A 

COSMOLOGICAL POPULATION OF TRIAXIAL 
HALOS 

2.1 Simulated cluster sample 

To begin, we generate a catalogue of triaxial dark matter halos with 
a generalised NFW density profile 



p{R) = 



R/R 3 (l + R/R 3 



(1) 



where S c is the characteristic overdensity of the halo, p c the critical 
density of the universe at the redshift z of the cluster, R 3 a scale 
radius, R a triaxial radius 



y-2 yZ y2. 

R2 = — + TT + — > (aO<c=l), 
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and a/c and b/c the minor:major and intermediate:major axis ra- 
tios, respectively. The NFW concentration is defined as usual 

R200 



C = 



Us 



and we employ a fully triaxial virial mass 

8007T 
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where ifeoo is the triaxial radius of the triaxial isodensity sur- 
face within which the mean density is 200 times the critical den- 
sity. We choose to define the virial mass as a fully triaxial quan- 
tity because such a definition is well-supported by collapse theory 
dSheth. Mo, & TormerJ hoOll) ) and is a more realistic representa- 
tion of dark matter halo shapes. 

We assemble a physically representative population of 15,000 
halos as follows 

• Draw a mass A/200 for each halo from the mass function of 
lEvrard et al. at z = 0.25, using the best fit parameters under 

the ACDM model for the differential number density of dark matter 
halos as a function of mass and redshift 



n(M,z) 



Ap m (z) 
M 



<Xeff(M) exp [-| ln<T _1 (M) 



■BY 



(5) 



where a e ff is the effective logarithmic slope and a 2 (M) is the 
variance of the density field smoothed on scales enclosing mass M 
at the mean density p m (z). We choose halos that range in mass 
from 1O 12,5 M to 1O 15,5 M in order to obtain a complete sample 
over the range of 10 13 — 10 15 Mq when potential scatter from 
lensing shape dispersion and triaxiality are accounted for. 

• Assign a concentration C to each halo according to the 
mass-concentratio n relationship for relaxed halos at z = 0.25 of 
iNeto et al] J2007I) . derived from the simulated groups and clusters 
in the Millennium simulation 



C = 



5.26 



M200 



W i h~ 1 M t 



with a log-normal scatter 



p(logC) = 



exp 



1 /logC-<logC) 



(6) 



(7) 



where (log C) is calculated via Eq.[6]for a given A/200 and the dis- 
persion a is taken to be 0.09 for masses less than 10 15 Mq and 0.06 
for mass es greater than tha t threshold, in rough agreement with the 
results of lNeto et al] d2007h . 

• Assign each halo a triaxial shape defined by minor and inter- 
mediate axis ratios a and b and two randomly distributed orienta- 
tion angles. The axis ratios are drawn fr om the distr i bution s found 
in the structure formation simulations of I Shaw et al. (2006). 

• Assign each halo a luminosity and richness drawn from the 
mass-richness (A^2oo : galaxy number) and mass-luminosity (L200 : 
i-band luminosit y) relations from the stacked weak lensing analysis 
of the SDSS of Ijohnston eTai] J2008l) : 



A^2oo (A/200) 



20 
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with 



A/ 



200 1 20 



■ ± QAstat ± 1-WO x 10 13 /i _1 M a 
a N = 1.28 ±0.04 



and 
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Figure 3. Offsets in the masses recovered, averaged across 50 noise realizations, for lensing by identical populations of 15,000 spherical and triaxial lensing 
halos, binned using SDSS observables. The black lines plot the percent offset of the lensing mass recovered from the triaxial lensing population compared to 
that for the spherical population in each luminosity L200 (left-hand panel) and richness N200 (right-hand panel) bin, for all 50 realizations. Overplotted are 
comparisons of the true mean values of the various mass estimators discussed in Sec. l2.2l to the true mean A/200 value for the triaxial population in each bin 
(all mass estimators are equal for the spherical population). Mes i s plotted with blue stars, M2DN FW with red crosses, and M^enfw with green circles. 
The effective spherical mass Mes is the triaxial mass measure most closely recovered by a stacked lensing analysis in low mass bins that contain many halos; 
the 2D lensing mass M2DNFW performs best in the high mass bins with few halos where shape averaging is incomplete. 



L 2 oo(M 200 ) = 40 ( W 10 h~ 2 L e (9) 

\™ 200 1 40 / 

with 

M 200 |4o = (9.5 ± 0A stat ± 1.2 sys ) x 1O 13 /i _1 M . 

a L = 1.22 ±0.04 

There is an intrinsic scatter in these mass-observable relations 
which is not well constrained but expected to be of order ~ 20%. 
In this study, because we seek to isolate the effects of triaxiality, we 
neglect this scatter and use the best fit values of the mass-richness 
and mass-luminosity parameters exactly with no error or scatter. 

For every case in this paper, an additional identical population 
of spherical halos with the same Af2oo, C, ^200. and L200 values 
is assembled for comparison. 

2.2 Measures of Mass 

Above, we defined M200 as a triaxial quantity, a choice we pre- 
fer because it describes halos more as we expect them to really 
be. However, many other mass measures have been used in previ- 
ous work, and the choice of measure can be important, as we will 
demonstrate later in this paper. Three alternatives for describing a 
triaxial lensing halo's mass are 

• An effective spherical mass and concentration where the ef- 
fective spherical virial radius r200.ES is defined as the radius of the 
spherical surface (which is not an isodensity surface) within which 
the mean density is 200 times the critical density. The effective 
spherical virial mass is the mass within that sphere 

M ES = (8007r/3)p c r? 00iES , (10) 

and the effective spherical concentration Ces is the ratio of 



r 2 oo,BS to the geometric mean of the triaxial scale radii Ces = 
f200, Es/r s ,ES, where r s , B s = R s (abc) 1/Z . Mes is always less 
than A/200, usually by only a few percent. This parameterisation is 
equivalent to directly calculating r2oo,BS from the raw mass distri- 
bution of a cluster in a structure formation simulation. This m ass 
estimator is very similar to that adopted by I Jing & Sutol J2002h to 
describe triaxial dark matter halo masses. Mes is better defined 
than Ces, because the concentration is meaningful only as a paral- 
lel to the NFW profile, where as the effective spherical mass is well 
defined for any model. 

• An effective lensing 2D circular NFW mass, calculated by fit- 
ting the projected density profile of a spherical NFW model to the 
projected density profile of the triaxial halo, giving mass and con- 
centration M2DNFW and C2DNFW ■ This mass makes sense only 
in a lensing context, as Midnfw will differ depending on the di- 
rection of projection for a given triaxial halo. 

• An effective 3D spherical NFW mass, calculated by fitting the 
3D density profile of a spherical NFW model to the density profile 
of the triaxial halo, giving mass and concentration Msdnfw and 
Czdnfw- This parameterisation is equivalent to fitting a spheri- 
cal NFW model to the mass distribution of a cluster in a structure 
formation simulation. 

All of these masses are equal for spherical NFW halos. 



2.3 Lensing Simulations 

Weak lensing distorts the shapes and number densities of back- 
ground galaxies. The shape and orientation of a background galaxy 
can be described by a complex ellipticity e s , with modulus \e s | = 
(1 — b/a)/(l + b/a), where b/a is the minor:major axis ratio, 
and a phase that is twice the position angle (f>, e s — \e s \e 2 "^. The 
galaxy's shape is distorted by the complex weak lensing reduced 
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shear, g — 7/(1 — k), where 7 is the complex lensing shear and 
k the convergence, such that the ellipticity of the lensed galaxy e 
becomes 



e +9 

1 



e +7 



(ID 



in the limit of weak deflections (the * denotes complex conjuga- 
tion). The distributions of ellipticities for the lensed and unlensed 
populations are related by 



d 2 e 



(12) 



assuming a zero-mean unlensed population, the expectation value 
for the lensed ellipticity on a piece of sky is (e) = g m 7. Thus 
the measured shapes of images can be used to estimate the shear 
profile generated by an astronomical lens. Lensing also changes 
the number counts of galaxies on the sky via competing effects; 
some faint sources in highly magnified regions are made brighter 
and pushed above the flux limit of the observation, but those same 
regions are stretched by the lensing across a l arger patch of s ky and 
so the number density of sources is reduced dCanizare i Jl982h ). 

The lensing convergence and shear of triaxial NFW halos are 
calculated numerically by integrating over combinations of the con- 
vergence of a spherical NFW and its derivatives, scaled by several 
factors that account for the distribution of mass along the axis ori- 
ented along the line-of-sight. A full de rivation of the lensing prop - 
erties of triaxial NFW halos is given bylOguri. Lee. & Sutol (12003), 
and extended in lCorless & Kind d2007l) . 

We simulate weak lensing through each of the 15,000 halos 
in each population, applying Eq.Qj]and the numerically calculated 
convergence and shear to lens catalogues of randomly positioned 
and oriented background galaxies with intrinsic shapes e s drawn 
from a Gaussian distribution with dispersion a = 0.2 in the mod- 
ulus je s j, number density of 5/arcmin 2 , and a redshift distribution 
modelled roughly on that of the galaxy population of SDSS Data 
Release 6 tha t peaks at z — 0.45 and has a tail out to z — 0.8 
dOvaizu et al.l {2007); note these are significantly lower redshifts 
than in a typical deep weak lensing study). For each triaxial popu- 
lation, the catalogue of background galaxies is randomly generated 
independently for every individual lensing halo. Background cata- 
logues exactly matching those of their corresponding triaxial halos 
are used for the halos in spherical population. Thus, we test the 
changes in the lensing signal and reconstruction when halos identi- 
cal in all but shape lens through an identical universe. 

All the lens halos are placed at the median redsh ift z = 0.25 
of the SDSS galaxy cluster population described in Koest eret al.l 
d2007l) . such that for the cluster population angular and physical 
scales become equivalent measures and the impact of the cluster 
redshift distribution need to not be accounted for. Again, this sim- 
plifying assumption is justified by our aim to isolate and understand 
the impacts of triaxiality. Magnification effects and Poisson noise in 
the galaxy distribution are accounted for. Throughout we assume a 
concordance cosmology with matter density parameter fi m = 0.3, 
Hubble parameter h — 0.72 in units of 100 km/s/Mpc, and a cos- 
mological constant Qa = 0.7. 



3 LENSING ANALYSIS 

We follow the binning procedure of Ijohnston et al.l d2008h and bin 
the halos in 12 richness and 16 luminosity bins. Within each bin 
we stack the lensed catalogues, assuming perfect knowledge of the 
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Figure 4. The residual mass offset in the recovered masses of the spheri- 
cal and triaxial halo samples for both luminosity (top panel) and richness 
(bottom panel) bins after correction by the effective spherical mass Me s ■ 
The triaxial bias is almost entirely removed in the low mass bins where 
the averaging over triaxial shape is most effective. The quantity plotted is 
{M 3ax - M S ph)/M S ph - ({M ES ) - <M 2 oo))/ (Afaoo). 



cluster centre, and fit spherical NFW models to the high signal- 
to-noise stacked data extracted from within an annulus covering 
0'.5 - 15'.0 (115 - 3500 kpc) from the field centre to estimate the 
mean mass of each bin. 

We primarily employ a radial shear profile fitting technique to 
fit the lensing masses. Under this method the tangential ellipticities 
et of the galaxies in the stacked catalogues are averaged in nine log- 
arithmic radial bins, giving a mean radial shear profile. This profile 
is then fit by minimising \ 2 



x 2 = 



£ 

i=l 



(et,< - g{ri-,U)) 2 



(13) 



where e t ,i is the mean tangential ellipticity in radial bin i, g is the 
reduced shear calculated for the spherical NFW model being fit at 
the mean radius of the ra dial bin rj ( f or the lensing properties of 
spherical NFW halos see lBartelman is the statistical 

error on e t ,i, and II is a two-element vector of the parameters defin- 
ing the model: virial mass M200 and concentration C. We assume 
a single effective mean lensing redshift of the background galaxies, 
which for our galaxy redshift distribution is z s = 0.435; for sim- 
plicity to isolate the effects of triaxiality we assume perfect knowl- 
edge of this value. It is possible to carry out a more precise analysis 
by scaling the measured shear of each galaxy via its redshift to a 
single source plane prior to radial averaging, but because we are 
not specifically interested in the impacts of the redshift distribution 
in this paper, but aim to isolate the effects of triaxiality, we adopt 
the simpler approach of assuming a single mean effective lensing 
redshift. The x is minimised over II using a simple downhill min- 
imiser from the GNU Scientific Library (GSL), and the values of C 
and M200 that minimise the function are returned as the mass and 
concentration estimate for the mass bin. 

To check that this method is robust (most especially our as- 
sumption that the lensed galaxies can be treated as if on a single 
plane at a single effective lensing redshift) we test our mass esti- 
mates using a maximum-likelihood method in which each galaxy 
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in the stacked catalogue is treated individually. We define the log- 
lik elihood function in the standard m ann er for weak tensing follow- 
ine lSchneider. King. & Erbenj d2000h and lKing & Schneiderl J200lh 

^ = -ln£ = -^lnp e (ei|g(^ ; n)) (14) 

i=l 

where u is the measured ellipticity of each galaxy, n 7 is the number 
of galaxies used, the reduced shear g is calculated for the spherical 
NFW model being fit, and n is again a two-element vector of the 
parameters M200 and C. We assign each lensed galaxy an individ- 
ual redshift, and again to isolate the effects of triaxiality, assume 
perfect knowledge of those galaxy redshifts in our analysis. The 
log-likelihood is minimised over II using the same GSL minimiser, 
and the maximum likelihood values for C and A/200 are returned 
as the mass and concentration for the mass bin. 

Normalizing against this more complex method we confirm 
that the averaged radial shear method provides accurate mass es- 
timates. Because the radial shear method is much faster computa- 
tionally than the maximum-likelihood method for large numbers of 
background galaxies (seconds vs. hours), we conduct the rest of our 
analysis using the averaged radial shear method only. 

These methods are different to those employed in 
lohnst onet alj |2008), which employs a non-parametric cross- 
correlation technique, but are very similar to those used in o ther 
weak lensing mass rec onstructions such as lClowe et alj J2006h and 
iLimousin et al.1 d2007h . We choose them for their simplicity and 
efficiency. Because we fit our NFW lens halos with NFW models, 
assuming a parametric model will not induce any error or scatter 
in our calculations because of disagreement between the profile 
shapes of the model and data. 



4 RESULTS AND DISCUSSION 

To begin, we simulate lensing by a single triaxial and correspond- 
ing spherical halo population (identical in all but shape) through 50 
different realizations of the universe: i.e. different randomly cho- 
sen background galaxy populations. To isolate the effects of tri- 
axiality, in each realization the triaxial and spherical populations 
lens through identical universes: the background galaxy catalogues 
are independently and randomly chosen for each lens halo in each 
realization, but are identical between triaxial and spherical coun- 
terparts. These 50 realizations of lensing through a single triax- 
ial and corresponding spherical population will effectively average 
over most of the intrinsic galaxy shape noise, and so we expect that 
for the spherical case the r ecovered lensing masses w ill accurately 
recover the true bin values. iMandelbaum et alj ( 120051) demonstrated 
that the masses recovered from a stacked lensing analysis in a uni- 
verse with a steeply sloped mass function, such as that simulated 
here, are expected to fall somewhere below the mean and above 
the median mass of the true bin population. If the mass dispersion 
within the bin is small, the difference between mean and median 
will be small, and that is indeed expected in this case in which we 
have included no intrinsic scatter in the mass-observable relations. 

Figure[T]plots the offset in the recovered bin masses, averaged 
across the 50 realizations, compared to the true mean and median 
mass values ((M200} and med(M2oo)) for the spherical and triax- 
ial halo populations, in both the luminosity and richness bins. The 
error bars show the standard deviation of the 50 realizations. 

To begin, we note that the masses of the spherical population 
(plotted with a solid black line) are recovered very well: indeed, 



the recovered masses tend to fall very slightly below the true mean 
values, and slightly above the true median values, exactly as ex- 
pected from the Mandelbaum result. The greater noisiness of the 
L200 masses is as expected; because there are four more L200 bins, 
the number of halos (and thus lensing signal) in each bin will de- 
crease and so the uncertainty on each individual realization will 
be higher. The higher levels of noise in the comparison to the me- 
dian mass values is also expected, because the median is a noisier 
quantity than the mean. Taken together, these indicate an accept- 
able recovery of the masses of the spherical population that follows 
the precedent of similar previous work; we are therefore confident 
to turn our attention to the recovery of the triaxial population. 



4.1 A triaxiality bias? 

The recovered masses for the triaxial population, averaged across 
the 50 realizations, are plotted with a blue dash-dot line in Figure 
[T] Crucially, the recovered mean mass values are consistently 1- 
2 percent less than those for their identical spherical counterparts. 
By this offset, which is found consistently in every realization, the 
lensing mass recovered from triaxial halos is shown to be a biased 
estimator of M2oo- The noise between the spherical and triaxial 
populations is correlated because we have lensed through identical 
universes with lens populations identical in every way except for 
their shape - thus the remnant galaxy shape noise is exactly the 
same for the spherical and triaxial cases. 

Though the triaxial bias is well within the statistical error bars 
in all bins, it is still problematic. Because the errors due to the shape 
noise of the lensing realizations are no longer centered on the true 
mass value - but are instead consistently offset by a few percent 
- comparisons between lensing constraints on the mass function 
and predictions from theory that assume symmetric or Gaussian er- 
rors will be incorrect. In particular, if the mean masses measured 
by lensing are consistently low, as we have demonstrated here that 
they are, then the mass function will be directly affected because 
the masses with which n(M) is calculated will be incorrect; this is 
crudely illustrated in Figure|2]in which the mean cumulative mass 
function N(> M) is plotted for the 50 realizations of Figure[T]for 
both spherical and triaxial populations. Note that this is indepen- 
dent of any change in the mass-observable relation used to sort ha- 
los into mass bins. As better surveys with more background sources 
and smaller measurement errors come online, this bias will become 
even more important. 

Why this apparent difference when triaxial halos lens through 
an identical universe? What physical or statistical explanation can 
be offered to explain why averaging does not work as expected in 
stacked lensing analyses? Figure [3] shows that the way in which 
halo mass is defined is crucial in understanding this apparent tri- 
axial bias. The black lines plot, for each of the 50 realizations, the 
percent offset of the lensing mass recovered from the triaxial pop- 
ulation M^ax compared to that recovered for the spherical popula- 
tion M Sp h 

M^ax — Ms p h n 
M Sph ' U3J 

in each luminosity L200 (left-hand panel) and richness A^oo (right- 
hand panel) bin. As seen before in the average in Figure [TJ the 
lensing masses consistently underestimate the triaxial masses M200 
compared to those of their matched spherical counterparts by ~ 
1 — 2% in the low mass bins and even more in the high mass bins. 
The very narrow dispersion between the 50 realizations illustrates 
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Figure 5. Offsets in the masses recovered, averaged across 50 realizations of the lensing populations, for lensing by identical populations of 15,000 spherical 
and triaxial lensing halos, binned using SDSS observables. The black lines plot the percent offset of the lensing mass recovered from the triaxial lensing 
population compared to that for the spherical population in each luminosity L200 (left-hand panel) and richness iVano (right-hand panel) bin, for all 50 
realizations. Overplotted in dashed red is the mean offset, and in solid blue is the mean value of the comparison of the true mean value of MeS to the true 
mean M200 value for the triaxial population in each bin, each averaged across all 50 population realizations. 
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Figure 6. Offsets in the masses recovered, averaged across 50 noise realizations, for lensing by identical populations of spherical and triaxial lensing halos in 
the 12th N20Q richness bin, plotted as a function of the number of halos in the bin, iV/, a ; . The left-hand panel plots the offset for each of the 50 realizations; 
the right-hand panel plots the dispersion around the mean deviation across the 50 realizations. The lower the dispersion, the more effective the shape averaging 
and the better a spherical stacked lensing mass estimate will recover the true bin mass. 



that the scatter in the triaxial bias due to galaxy shape noise is far 
less than the bias itself, indicating the bias is robust in this respect. 

Overplotted in the figure are comparisons of the true mean al- 
ternate masses - calculated via the alternate mass measures listed 
in Section l2^2l - to the true mean triaxial mass for the triaxial popu- 
lation in each bin (the mass estimators are all equal for the spheri- 
cal population): (M a it — M2no)/M2oo. It is hoped that the lensing 
mass recovered from the triaxial population is a good estimator of 
one of these alternate mass measures. And indeed, the recovered 
lensing masses are seen to better estimate all three alternate mass 



measures. Of the three, the 3D NFW mass Msdnfw, plotted in 
green, provides the least improvement, as the lensing mass consis- 
tently overestimates it by a small percentage. 

However, in the low mass bins the recovered lensing mass ap- 
pear to be an excellent estimator of the effective spherical mass 
Mes (Eq. I lOt . and in the higher mass bins of the 2D NFW lens- 
ing mass M2dnfw- In those higher mass bins where M^dnfw 
is best recovered there are only a few halos (< 100) to average 
over; in these cases the 2D lensing mass estimator performs best 
because, though it still demands averaging over ellipticities in 2D 
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to circular symmetry, it, by its 2-dimensional definition, does not 
require averaging over the unseen third dimension along the line 
of sight. However, such a 2D mass estimator is of less general ap- 
plicability because it is meaningful only in the context of lensing 
and requires the adoption of a line-of-sight orientation to be calcu- 
lated. Theory makes predictions for 3D mass functions - and it is 
only through comparison to these predictions that a measured mass 
function can be leveraged to constrain cosmological parameters. 
Translating predictions to a derivative 2D mass is undesirable be- 
cause their strength is diminished, though in problems in which the 
accurate estimate of the highest mass bins is critical those disad- 
vantages may be outweighed by its accuracy. When at all possible, 
accurate 3D masses are preferable. 

Luckily, in the lower mass bins where the shape averaging 
is effective, the lensing mass performs even better as an estima- 
tor of the 3D effective spherical mass. This is further illustrated 
in Figure Q] where the lensing masses recovered from the triaxial 
population, averaged across the 50 realizations, are compared to 
the true effective spherical mean and median mass values ({Mes} 
and med(MEs))- In the low mass (many halo) bins, the errors in 
the stacked lensing masses recovered for the triaxial population ex- 
actly match those for the spherical population. This is confirmed in 
Figure [4] which shows the residual offset between the triaxial and 
spherical lensing masses when the true Mes values of the triaxial 
halo population are used to remove the triaxial bias; they are very 
close indeed to zero in the low mass bins where averaging is ef- 
fective. We will return to the important question of when averaging 
over triaxial halo shape is effective in Section l4~2l 

Now, to test the robustness of the apparent triaxial bias across 
different realizations of the lensing halo population, we perform 
lensing through 50 different halo populations, each chosen ran- 
domly from the distributions outlined in Sec. 12. II Again, identical 
triaxial and spherical populations are generated for each realization, 
and the background lensing catalogues are independently and ran- 
domly chosen for each lens, though matched between triaxial and 
spherical counterparts. Figure[5]again plots the percent offset of the 
lensing mass recovered from the triaxial population compared to 
that for the counterpart spherical population for all 50 realizations, 
together with the comparisons of the alternate mass measures to 
M200, as plotted in Figure [3] The triaxial bias is present as before, 
with greater scatter due to the noise introduced by the 50 different 
halo populations. The scatter is very large in the high mass (few 
halo) bins where the averaging over triaxial shape is insufficient, 
because in these bins the shape distribution is quite different for 
every realization. The scatter in these bins is symmetric, showing 
that the large underestimation of the mass seen in Figure ?? for the 
single lens population is only one possible outcome when the num- 
ber of halos is too small; insufficient averaging can lead to both 
under- and over- estimates of the effective spherical mass. Again, 
Mes, whose average across the 50 realizations is plotted in blue 
in Figure [5] is the mass measure best estimated by the recovered 
lensing mass, as long as the averaging is sufficient. Our findings 
are therefore shown to be robust both for different realizations of 
the lensing universe and of the lensing population. 

Though previous work in ICorless & Bngl d2008l) has shown 
averaging over spherical models individually fit to weak lensing 
data to give approximately correct mean values for triaxial mass 
M200 and concentration over a population of 100 equal mass triax- 
ial halos, it was not concerned with errors of a few percent. Those 
finding are therefore consistent with those in this paper, in which 
the averaging problem is examined in far greater detail. 

Though the bias seen in the recovery of the triaxial mass func- 



tion is small, at this time significantly smaller than other sources 
of error (crucially in the normalisation and scatter in the mass- 
observable relationships), it is important. In future surveys such as 
those reviewed by the DETF its relative weight compared to other 
errors will increase as those other sources of error are controlled. 
Even now, a bias of any size that consistently underestimates mass 
can be important given the exponential dependence of the normal- 
ization of the m atter power spectrum on the high-mass end of the 
mass function jAlbrecht et alj|200^) ; given that the bias is easily 
removed, it should be. 

We have demonstrated that the presence of triaxial halos in a 
lensing population can affect the measurement of average masses in 
bins typical of those used in stacked analyses to measure the mass 
function. The choice of mass becomes important, as the presence 
of triaxiality breaks the degeneracy between many possible mass 
definitions. We showed that the 2D lensing mass is best recovered 
when shape averaging is not adequate, but that the 3D effective 
spherical mass is best when there are enough halos in a given bin 
to make the averaging over triaxial shape effective. 



4.2 Effective Triaxial Averaging 

But how many halos are required for that to be so? To test this, we 
take a single mass bin - the 12th richness bin, at the uppermost end 
of the mass function - and populate it with an increasing number 
of halos N. We do this for 50 realizations of paired triaxial and 
spherical halo populations. We are looking for the number of halos 
at which the 'triaxial bias' becomes constant - i.e. the relationship 
between the recovered masses for the triaxial and spherical popu- 
lations is constant, such that it can easily be accounted for by using 
Mes- The left-hand panel of Figure [6] plots the percent deviation 
between the triaxial and spherical recovered masses as a function of 
for all realizations, and the right-hand panel plots the dispersion 
of those deviations around the mean over the 50 realizations. As ex- 
pected, the dispersion is largest when the number of halos in the bin 
is smallest. It decreases significantly to about 0.3% by N=500, and 
flattens out to 0.2% by N — 1600. This plateau value represents 
the dispersion due to a combination of the noise of the different re- 
alizations of the lensing population and the intrinsic galaxy shape 
noise. 

When designing strategies for present and future surveys, it is 
important to ensure there are an adequate number of halos in each 
mass bin for averaging over triaxial shape. Our results suggest 100 
halos as a minimum number, and 500 as highly preferable. For the 
upper reaches of the mass function where such numbers will not be 
available, i ndividual triaxial mass m odelling using techniques such 
as those in ICorless & Kingl J2008I) may be necessary, as the ap- 
proximation of average sphericity no longer holds. Neglecting this 
caution may result in errors of up to ~ 5% in virial mass estimates. 



4.3 Normalization and Scatter in the Mass-Observable 
Relations 

In addition to introducing a potential bias in direct estimates of 
masses from statistical weak lensing, the presence of triaxial struc- 
ture in dark matter halos also can introduce an effective scatter into 
the mass-observable relations that are used to sort halos into bins 
for stacked analysis. To illustrate this, Figure|7]plots the bin mem- 
bers of the 6th, 7th, and 8th richness bins for a single triaxial halo 
population. The top panel plots them as function of M200, the mid- 
dle as a function of Mes, and the bottom panel as a function of 
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Figure 7. Distributions of triaxial masses M200 (top panel), 3D effective 
spherical masses Mes (middle panel), and 2D circular lensing masses 
M2DNFW (bottom panel) of the halos in the 6th-8th richness bins for a 
single triaxial population. The standard deviation within each bin increases 
by 1 — 5% if the lensing mass M2DNFW i s used rather than the triaxial 
or effective spherical masses. 

M2dnfw- For a spherical population all of these would be identi- 
cal to the M200 distribution (note that we neglect all intrinsic mass- 
observable scatter). 

In this simulation we have assumed that the true mass- 
observable relation is between the observable and the triaxial halo 
mass M200 when sorting halos into mass bins. This may not be 
the case; perhaps the tightest mass-observable relation with least 
scatter utilizes the effective spherical mass, or some other mass 
measure. However, what this figure clearly illustrates is the critical 
role consistent and well-understood mass estimators play in accu- 
rately characterizing the mass-observable relations. Quantitatively, 
the scatter om in the mass-richness and mass-luminosity relations 
using the the mean effective lensing mass (Midnfw) varies from 
about 6% in low mass bins to over 8% in high mass bins, compared 
to no scatter in the triaxial mass Af2oo- For Mes, the effective scat- 
ter is ~ 2%. 

As an example, determining the scatter in the mass-observable 
relations by modelling individual galaxy clusters using a standard 
spherical lensing model, which gives an estimate of M^dnfw, 
will almost certainly significantly increase the effective scatter in 
the mass-observable relation, this is because there is little reason to 
expect richness or luminosity (or X-ray luminosity or temperature) 
to be strongly dependent on lensing orientation, and so comparing 
these quantities with 2D lensing masses will necessarily increase 
the scatter in the apparent mass-obse rvable relations. Such scat- 
ter is widely known to be important jAlbrecht et all J2006h : Ivoil 
liool)) because the mass function is steeply sloped, such that there 



are more low mass halos to scatter up than high mass halos to scat- 
ter down; this causes the mass recovered for a bin to decrease as 
the scatter increases. The scatter among bins must be well known 
in order to correct for this effect and accurately interpret the mass 
measurement in each bin. 

Regarding the mass-observable normalizations, if a stacked 
lensing analysis such as this is used to determine the normaliza- 
tion of the mass-observable relation, it is crucial to recognize that 
it is the relationships of observables to Mes that are being char- 
acterized, and not those to M200 or any other mass estimator. If 
care is not taken to identify and account for the 3D mass estimator 
that the lensing analysis accurately recovers, future studies that use 
the mass-observable relations derived from lensing to constrain the 
mass function directly from observables without any independent 
measure of mass will suffer from the same mass biases of the weak 
lensing analysis. 

At present, the increase in mass-observable scatter due to tri- 
axiality, even using the poor M2DNFW indicator, is small com- 
pared to the best estimates of the intrinsic scatter in most mass- 
observable relations. Recent work puts the intrinsic scatter in the X- 
ray mass-temperatu re and mass-lum i nosity r elations at ~ 1 5% an d 
45% respectively jVikhlinin etail J2006I) ; IStanek et all {2006)), 
and that in the mass-richness relation perhaps as h igh (the relation 
betwe en richness and X-ray luminosity is found bv lKochanek et al] 
(2003) to have scatter of 0.33 dex, greater than 100% for some 
values of L x ). Most optimistically, the relation between M500, a 
mass measure frequently used in X-ray analyses, and Yx , a proxy 
obtained by multipl ying together the gas mass and X-ray temper- 
ature was shown in iKravtsov. Vikhlinin. & Nagaa ( feOOot) to have 
an intrinsic scatter of less than 10%. Thus, in present analyses us- 
ing all but the last relation - which is unavailable as an observable 
for large optical weak lensing surveys unless they are accompanied 
by extensive complementary X-ray observations - the scatter due 
to triaxiality is non-negligible but unlikely to be of primary con- 
cern. However, in future surveys in which the intrinsic scatter in 
the mass-observable relations can be controlled to a few percent, as 
the DETF argues it must be, then the scatter due to triaxiality must 
be carefully considered and controlled. 

4.4 Summary 

We investigate the impacts of triaxiality on the accuracy of mass es- 
timates derived from stacked weak lensing data with the aim of con- 
straining the cluster mass function. Though the quantitative results 
of our study may not directly translate to other analysis methods 
using non-parametric reconstructions or cross-correlation stacking 
techniques, the general finding will hold true: the presence of triax- 
iality makes the choice of mass estimator crucially important. 

Most importantly, we find that triaxiality can introduce a bias 
in the masses recovered from a triaxial halo population using a 
stacked lensing analysis, giving mass values that are a few percent 
low in all bins compared to those of their spherical counterparts. 
This bias in the mean masses statistically measured by lensing will 
lead directly to incorrect constraints on the mass function n(M). 
Further, calibrations of the scatter in, and the normalization of, the 
mass-observable relations critical to accurate measurement of the 
mass function derived from lensing will also be made inaccurate by 
the bias. However, we show that this bias is removed if the effective 
spherical mass Mes, rather than the triaxial virial mass M200, is 
employed in defining the triaxial halo masses. While the effective 
spherical mass M es is not the ideal mass measure for all problems 
- especially those studying individual clusters (see lCorless & Kingl 
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(2008))- and has the disadvantage of treating dark matter halos 
as approximately spherical when all predictions suggest a signfi- 
cantly triaxial population, it is most suitable in this case. Further, 
we demonstrate that least 100, and preferably 500, halos are needed 
in every mass bin to make averaging over triaxial shape effective 
and so reliably recover the effective spherical mass Mes- 

Thus, stacked lensing analyses aiming to constrain the mass 
function or characterize the mass-observable relations provide ac- 
curate mass estimates only if the 3D mass estimators employed are 
very carefully selected and an adequate number of halos are stacked 
in each observable bin. With these caveats, applied by carefully cal- 
ibrating any stacked lensing mass reconstruction method with sim- 
ulations, they can be applied to future weak lensing surveys such 
as JDEM, LSST, and DES unhindered by galaxy cluster and group 
triaxiality. 
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